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Abstract 

A natural formulation of the theory of quantum measurements in 
continuous time is based on quantum stochastic differential equations 
(Hudson-Parthasarathy equations). However, such a theory was developed 
only in the case of Hudson-Parthasarathy equations with bounded coef- 
ficients. By using some results on Hudson-Parthasarathy equations with 
unbounded coefficients, we are able to extend the theory of quantum con- 
tinuous measurements to cases in which unbounded operators on the sys- 
tem space are involved. A significant example of a quantum optical system 
(the degenerate parametric oscillator) is shown to fulfill the hypotheses 
introduced in the general theory. 



1 Introduction 

A powerful formulation of the quantum theory of measurements in continuous 
time is based on quantum stochastic calculus |20I21| . In such an approach, 
the quantum stochastic Schrodinger equation, or Hudson-Parthasarathy equa- 
tion (HP-equation), is combined with suitable field observables [31517^ 9]: the 
resulting formulation is particularly suited for applications in quantum optics 
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and for building up a photon detection theory |4I5I11I18I25| . However, the the- 
ory is fully developed only for the case in which the HP-equation involves only 
bounded operators in the initial Hilbert space. Many results are known on the 
existence and uniqueness of the solution of the HP-equation with unbounded 
coefficients [13115116] : our aim is to combine these results with the equations 
for the (unbounded) field observables and to show how to arrive to the key 
evolution equation (j40p of the theory of continuous measurements, which con- 
cerns the "reduced characteristic operator" (P7)) . Moreover, for applications, 
it is important to consider the case in which the initial state of the quantum 
fields is not only the vacuum, but at least a generic coherent vector. This gives 
that the reduced dynamics is not a quantum dynamical semigroup and we need 
to handle a master equation with a time-dependent, unbounded Liouville oper- 
ator. Finally, a relevant physical example is given: the degenerate parametric 
oscillator [11]. The paper is based on Castro's PhD thesis [12] . 

For any separable complex Hilbert space f) let us introduce the following 
classes of operators on it: £(()), the space of bounded linear operators, U{t)) the 
class of the unitary operators, T(f)) the trace-class, §([)) := {p G T(f)) : p > 0, 
Tr{p} = 1} the set of statistical operators. 

Then, we introduce the symmetric Fock space 3" over L^(R_|_;Z), where Z 
is a d-dimensional complex Hilbert space (the multiplicity space) in which we 
fix a complete orthonormal system {zi, i = l,...,ci}. We denote by e(/) 
the exponential vector in the Fock space associated with the test function 
/ g L^(IR_|_;Z) and we call coherent vector ■ip{f) := ||e(/)||~"^ e(/). Recall that 
{e{g)\e{f)) = exp(g|/). We assume familiarity with such notions and with 
quantum stochastic calculus [H]. We shall use the notation fi{t) := {zi\f{t)) 
for all z > 1 and we set fo(t) = 1. We fix the sets 

M = L^{R+;Z)nL^^{R+;Z) and linear span of {e(/) :/ e M}. 

The set £ is dense in J Qli Corollary 19.5 p. 127]. 

An important feature of the Fock space 3^ is its structure of continuous 
tensor product. For any choice of the times < s < t let us introduce the 
symmetric Fock space 3^(s,t) over L'^{{s, t); Z) and the symmetric Fock space 3^(t 
over L'^{{t, oo); Z). Then, we have the natural identifications 

J~ 5-(o,3) «) J(,,t) ® J(t and e(/)~e(/(o,,))®e(/(,,t))(g)e(/(t), (1) 

where f(s,t){x) ■= l(s,t){x)f{x) and /(^(x) := l(t,oo)(a;)/(a;). The symbol (g) de- 
notes the tensor product of Hilbert spaces, vectors and operators; the algebraic 
tensor product of dense spaces is denoted by 0. 

The Weyl operator j21J W{g; U), with g G L^{m.+ -Z) and U e U{L'^{R+; Z)), 
is the unique unitary operator on 3^ defined by 

Wig;U)e{f) = exp[-^\\gf^{g\Uf)}e{Uf + g), yfeL\R+;Z). (2) 
From the definition one obtains the relations 

W{g; Uy' = W{g; U)* = W{-U*g- U*) 
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and the composition law 

W{h; V)W{g; U) = exp { - ilm{h\V g)]W {h + Vg; VU). (3) 

Finally, we denote by Ajj(t), Aj{t) the creation, gauge and annihila- 

tion processes associated with the basis {z^, i = 1, . . . , d}; we shall use also the 
notation 

A,o(i) = 4W, Ao,(t) = Aoo(i)=^, i,j = l,...,d. (4) 

In particular we have (e(g)|Ay (t)e(/)) = 9i{s) f j{s)ds {e{g)\e{f)) , i,j = 
0,...,d. 

Let 3-C be a separable complex Hilbert space, the initial space, and let us call 
S(K the quantum system with Hilbert space 

We refer to [15|21j for the definition of quantum stochastic integrals with 
respect to the operator noises A^ , but we need to report at least the notions of 
adapted process and stochastic integrability. 

Definition 1 ([3TJ p. 180], [1^1 Definition 2.1]). Let _D be a dense manifold in 
J{. A family {L{t),t > 0} of operators in 3-C (g) is an adapted process with 
respect to (£>, M) if (i) D £ C nt>o Dom(L(t)), (ii) the map t H> L{t)u (g) e(/) 
is strongly measurable, G D, / e M, (iii) L{t)u e(/(o,t)) G (g) 3^(0.*) and 
L{t)u e(/) = {L{t)u e(/(o,t))) ® e(/(0, Vt > 0, u e I?, / e M. 

If additionally the map 1 1— >■ L(t)u ® e(/) is continuous for every u € D and 
/ e M the process is said to be regular adapted. Moreover, the adapted process 
L is said to be stochastically integrable if, for alH > 0, m G D and / G M, one 
has ||L(s)u (g) e(/)|pds < cxj. 

A key notion in the construction of dilations of quantum dynamical semig- 
roups is the one of cocycle [1]. We introduce the strongly continuous one- 
parameter semigroup {d{t), t > 0} of the shift operators on L^(M+; Z) and its 
second quantisation Q on S': for every t > 

{9tf){x)^f{x + t) and Q,e{f) ^ e{9tf), yfeL\R+;Z). (5) 

Let us note that, for r < s, (^9tl(^r,s)) i^) — ^{r,s){x + t) — l(^r-t,s-t){x)i this 
implies 

0te(/(o,s)) = e(0), for < s < i. Of C S'^r-t^s-t) for < t < r < s. 

Moreover, it turns out that 6^* is an isometry. We extend 9* to the space JC®? 
by stipulating that it acts as the identity on 

Definition 2 (Right and left cocycles). A bounded, adapted operator process 
X(t) in 3-C(g)S' is called right cocycle (respectively, left cocycle) if for every s,t > 
we have X{t + s) = e,*X(t)e,X(s) {X{t + s) = X {s)e * X {t)e . 
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2 The Hudson-Parthasarathy equations 

Let us consider the quantum stochastic differential equation (QSDE) for oper- 
ators on 5f (g) 3^, known as right HP-equation: U{0) = 1, 

d'^(^) = ( J2 R^^Mit) + J2 ^*jdA'j (i) + J2 + Kdt^ U{t), (6) 

i>l i,j>l i>l 

where the coefhcients K, Ri, Ni, Fij, with i,j = l,...,d, are (possibly un- 
bounded) operators in the initial space IK. Very general sufficient conditions, 
which guarantee the existence of a unique solution of (jH]) and the fact that such 
a solution is a unitary cocycle, are given by Fagnola and Wills |16| . 
By using the notation ^ and by setting 

Fao = K, F,o = R, , Foj = Nj , (7) 

we can write the right HP-equation in the shortened form 

dUit) = J2 (t) Uit), C/(0) = i. (8) 

We shall need also the adjoint equation, the left HP-equation: 

dV{t) = V{t) Fn ^(0) = 1. (9) 

i,3>0 

Definition 3 (Right Solution - [151 Definition 3.2]). Let I? be a dense subspace 
in 5{. An operator process C/ is a solution of the right HP-equation in 13 £ for 
the matrix F if: 

(i) each operator 11 is closable and |J U{t){D £) C Q Dom(Fij 1); 

t>0 ij^O 

(ii) each process Fij 1 [/ is stochastically integrable and 

U{t) = \+Y^ [ 1 C/(s)dAy(s) onD0£, Vt > 0. 

Definition 4 (Left Solution - [15l Definition 3.1]). Let D be a dense subspace 
in 5{. An operator process is a solution of the left HP-equation in 13 £ for 
the matrix F* if: 

(i) D C ni,j>oDom(Fy*) and the linear manifold (Uij>o^y* (-^)) © is 

contained in the domain of V(t), Vt > 0; 

(ii) the processes (y{t)F^*;t > O) are stochastically integrable and 

V{t) = 1 + 5Z / ^i^)Fji dAij(s) on I) £, Vi > 0. 
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Hypothesis 1. (The matrix F) 

(i) F = {Fij;0 < i,j < d) is a matrix of closed operators in the initial space 

^K. By F* we denote the adjoint matrix, defined by {F*).ij — F*^. We also 
define Dom(F) := n,,,>o Dom(^^,,), Dom(F*) := n.,,>o Dom(F,* ) 

(ii) For I < i, j < d, we have Fij = Sij—Sijl, where the Sij are bounded operat- 

ors on D-C satisfying the unitarity conditions X]fe=i ^ki^kj — Y^t=i ^ikS*/. = 
Sij. 

(iii) There exist a dense subspace D which is a core ioi K, Ri, Ni, i — 1, . . . ,d, 
and a dense subspace D which is a core for K*, R*, N* , i = 1, . . . , d. 

(iv) Dom(7V*) D DUD, Dom(i?i) D DUD, Dom(7Vi) D T)om{K), \/i > 1. 

(v) Vfc > f , Vm e Dom(X): SkiU e Dom(i?*), > 1. 

(vi) The operators K and K* are the infinitesimal generators of two strongly 
continuous contraction semigroups on J{. Moreover, we have Vu G D, 
\fveD 

2Re{Ku\u) = -^\\Rku\\^, 2Re{K*v\v) = -J2\\Nkv\\^ ■ (10) 

k>l k>l 

(vii) N*u = - J2k>i S*MRkU , Wu€DuD, Wi>l. 

(viii) There exist a positive self-adjoint operator C on 5£ and the constants 
(5 > and 6i, 62 > such that [H pp. 281-291] Dom{C^^^) C Dom(i^) 
and 

(a) for each e e (0,(5), there exists a dense subspace D^ C D such that 

Cy^L», C D and each operator F*jCl^^\ is bounded, where = 
C 

(b) for all < e < (5 and uq,. . . ,Ud € Dom(F), the following inequality 
holds: 

{Ui\CeFijUj) + {FjiUi\CeUj) + '^{FktUi\CeFkjUj) 
ij>0 \ k>l J 

<Y,{bl{u^\CeU^)+b2\\u,\\^) . 
i>0 

Proposition 1. Under Hypothesis{Ii also the following properties hold: 

1. Dom(i?fc) D Dom(A') U Dom(i^*), Dom(iV*) D Dom(X) U T)oia{K*), 
k>l. 

2. Dom(i^) = Dom(ii:); Eqs. ^ holdWu € Dom(i^), Vw € Dom{K*). 
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3. N*u = - ^^>i Sl^RkU, R,u = - Efe>i S^uN^u, Vu G Dom(i^)UDom(i^*), 
Vi > 1. 

I N,u ^ - J2k>i RlSk^u , Vw G Doni(i^), = 1, . . . , d. 

5. /or every choice of uq, ui, . . . , Ud in Dom(i^) and of vo, wi, . . . , Vd in D , 
we have 



{u,\F,jUj) + {Fj,u^\uj) +J2{Fk^u,\FkjUj) j = 0. (11a) 

i,j>0 \ k>l 

E + (^.».> + E iP^*k^^\P,>l) 1 ^ 0- (lib) 

i,j>0 \ k>l 

Proof. By ^ we get, V0 G D, \\Rk(p\\'^ < 2\{K(p\(j))\. For any u G Doni(i^) we 
can find a sequence u„ G -D converging to u. Then, Kun — >■ i^it weakly and by 
the proposition at p. 112 of 23 the sequence Kun is norm bonded: ||i4rw„|| < c. 
Then, we have \\Rk{un - Um)|P < 2|(if(u„ - Mm)|u„ - Um)\ < 4c||w„ - UmW- 
Therefore, i?fcU„ is a Cauchy sequence. Being Rk closed, {u, lim„ RkUn) belongs 
to the graph of Rk and, so, u G Doni(i?j.); this gives Doni(i?j.) D Dom(_fir). By 
the same property of closed operators, we get that the first equation in (jlOp can 
be extended to the whole Dom{K). Similarly, from the second equation in (|10l) 
we get that it can be extended to the whole Dom{K*) and that Dom(A^^) D 
Dom(i<'*). Once again by the property above of closed operators and by the 
unitarity of the operator matrix S, we get that point (vii) can be extended 
to Dom{K) U Dom(_ftr*). Therefore, on the same domain, J2k>i ll-^fc""!!^ — 
J2k>i ll-^fc^ll'^- exchanging R^ and in the two equations in pUj) we get 
also Dom(i?fc) D Dom(ii'*) and Dom(A^^) D Dom(_ft'). By these results and the 
definition of Doni(i^) we get Dom(i^) = Dom(i^). This ends the proof of points 
(l)-(3). ^ 

By using the first equation in point (3) of this proposition and point (v) in 
Hypothesis m we have, Vw G Dom{K), Vu G Dom(ii') U Dom(7^*), 

{V\N,U) = {N*V\U) =-Y,{S*krRkv\u) ^-J2{Rkv\Sk^u) = -J^ {^iKSkru) . 

k>l k>l k>l 

By the density of Dom{K) U Dom(/\ *) in J{ we have the statement in point (4). 
Equations in point (6) are by direct verification. □ 

Note the difference in the domains of Eqs. (|lip . The last requirement in 
point (iv) has been added just to have Eq. (|llal) on the whole Dom(F) and not 
only on D. As one can check, also Hypothesis HGC at p. 205 of [15^ holds for 
the coefRcients of the left HP-equation ([9]) . We collect in the following theorem 
many results. 



6 



Theorem 2 ([HI Proposition 2.2, Theorem 2.3]; [TH Proposition 6.3, The- 
orems 8.4, 8.5]). Under Hypothesis [7] the left HP-equation ([9]) has a unique 
solution {V{t)] t > Q\ on D Q E, which is a strongly continuous left cocycle of 
contractions. Moreover, U{t) — V(t)* is a right cocycle and solves the right 
HP-equation on Dom(C^/^) £. 

As we take U = V* , ii V is an isometry process, C/ is a coisometry process 
and vice versa. The fohowing Proposition is a small variation of Corollary 2.4 
of [H] or of Corollary 11.2 of [H]. 

Proposition 3. Under Hypothesis{^the contractive solution U of ^ introduced 
in Theorem \^ is a strongly continuous isometric process. Moreover, if U is 
unitary, it is the unique bounded solution on Dom(C^^'^) £ of ([8]). 

Proof. Let U be another bounded solution and apply the second fundamental 
formula of QSC ([21 , Proposition 25.2], ^ Eq. (2)]) to U, U. We get, V/,5 G M, 
Vu,w e Dom(Ci/2), 

{U{t)v e{g)\U{t)u e(/)) - (w e{g)\u(E) e(/)) 

= V / ds g,{s){{U{s)v(g)e{g)\Fij(g)tU{s)u(g>e{f)) 

+ {Fj, ®lU{s)v®e{g)\U{s)u®e{f)) 

+ {Fk^ 1 Uis)v e{g)\Fkj U U{s)u e(/))}/j(s) . 
fe>i 

But U and U are solutions on Dom(C^/^) £; by point (i) of Definition [3] 
and Eq. (|lla|) . we have that the integrand vanishes. Therefore, V/, g G M, 
\/u,v € Dom(Ci/2), {Uit)v (g) e{g)\U{t)u(g) e{f)) = {v (g) e{g)\u (g) e{f)) and 
this is equivalent to U{t)*U{t) = 1. This equation ioT U — U gives that U{t) 
is isometric, while for U unitary gives U{t)* = U{t)*. Being the adjoint of a 
strongly continuous process, U is weakly continuous and, being an isometry, it 
is strongly continuous. □ 

3 Quantum dynamical semigroups and unitary 
cocycles 

For the physical interpretation of the evolution operator U{t) we need it to 
be a strongly continuous cocycle of unitary operators (see [5] Section 2.2 and 
references there in). The unitarity is associated to some property of a related 
quantum dynamical semigroup (QDS); so, we start with some notions on QDSs. 

Definition 5 (QDS). Let us consider a family {T{t), t > 0} of bounded oper- 
ators on £(IH) with the following properties: 
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(i) T{t + s)^ T{t)T{s), ys,t> 0, and T(0) is the identity map; 

(ii) T{t) is completely positive, Vi > 0; 

(iii) T(t) is a cr-weakly continuous operator on L{'K)^ Vt > 0; 

(iv) for each X e L{JC) the map t !—> T(t)[X] is continuous with respect to the 
cr-weak topology of L{JC). 

Then, the family of operators T{t) is called a quantum dynamical semigroup. If 
also T(t)[l] = 1 holds Vi > 0, the QDS T{t) is said to be Markov or conservat- 
ive. 

Theorem 4 ([Ml Theorem 3.22, Corollary 3.23]). Let A be the infinitesimal 
generator of a strongly continuous contraction semigroup in 'K and let L^, k = 
1, . . ., he operators in such that the domain of each operator contains the 
domain of A and for every u G Dom(y4) we have 2Tie{u\Au)+^f,y^ ||Lfeu||^ = 0. 

For all X g 'C(5f), let us consider the quadratic form C[X] in 'K with domain 
Dom(^) X Dom(^) given by 

{v\C[X]u) = {v\XAu) + {Av\Xu) + ^{Lkv\XLku). (12) 

k>l 

Then, there exists a QDS T{t) solving the equation 

{v\T{t)[X]u)^{v\Xu)+ f {v\C[T{s)[X]\u)ds (13) 

Jo 

with the property that T{t)[l.] < 1, Vt > 0, and such that for every a -weakly 
continuous family T'(t) of positive maps on £(Jf) satisfying Eqs. (jl2[) and (jl3p 
we have Tit)[X] < V{t)[X], Vt > 0, for all positive X e £(Jf). 

// moreover the QDS T{t) is conservative, then it is the unique a-weakly 
continuous family of positive maps on £(Jf) satisfying Eq. il3\} . 

The QDS T(t) defined in Theorem |4] is called the minimal quantum dynam- 
ical semigroup generated by A and Lk, k = 1, . . .. Sufficient conditions to assure 
Markovianity of a QDS are known [13ll5j . In the application we shall use the 
following result. 

Theorems ( |15[ Theorem 9.6]). Let A, Lk be as in Theorem^and suppose that 
there exist two positive self-adjoint operators Q and Z in 'K with the following 
properties: 

• T)om{A) is contained in T)QTn{Q^) and is a core for ; 

• the linear manifold p|j,>-^ (Dom(^^)) is contained in I}om{Qi) ; 

• Dom(A) C Dom(Z5) and 

-2Re{u\Au) = \\Lku\\^ = \\zK\\^, Vu e Dom(A); 

k>l 
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• Dom((5) C Dom{Z) and for all u G Dom((52) we have \\Z'^u\\ < \\Q2u\\; 

• there is a positive constant b depending only on A, Lk, Q such that, for 
all u G Dom(A^) , the following inequality holds 

2Re{Q^u\QiAu) + ^ WQ^kuf < 6|lQ^uf . 

A:>1 

Then, the minimal quantum dynamical semigroup associated to A and Lk is 
Markov. 

Now, we can go back to the problem of the unitarity of U{t). 

Theorem 6 ( [151 Theorems 10.2, 10.3]). Under Hypothesis [I\ the contractive 
left cocycle V solving ([9]) is such that the family of operators T(t) defined by 

{v\f{t)[X]u) = {V{t)v®e{0)\{X(g,l)V{t)u(g)e{0)), Vu,wGM, yX e L{3i). 

is the minimal QDS generated by K* and N^, k — I, . . . , d. 
Moreover, the following conditions are equivalent: 

(ii) the process V is an isometry; 

(ii) the minimal QDS associated with K* and is conservative; 

Hypothesis 2 (Markov condition). The minimal QDS generated by K* and 
N^ , k — 1, . . . , d, is conservative. 

Now, both U{t) and V{t) ~ U{t)* are isometrics and, so, they are unitary 
operators. 

Corollary 7. Under Hypotheses]^ and\^ the process U introduced in Theorem 
is unitary and it is the unique bounded solution on Dom(C^/^) £ of ([5]). 

Remark 1. When U{t) is a strongly continuous unitary right cocycle, one can 
define the unitary evolution 

U{t, s) := U{t)U{s)*, t>s>0. (14) 

It is easy to check that it is strongly continuous in t and s and such that 

U(t,s) = Q;U{t-s)es, U{t,r) ^U{t,s)U{s,r), 0< r<s<i. (15) 

Moreover, the operator U (t, s) is adapted to Jf ?(s.t) in the sense that it acts 
as the identity on 9^(o,s) ® 3^(t and leaves Jf ® 3^{s,t) invariant. The unitary 
operator U{t, s) is interpreted as the evolution operator of system Sj{ and fields 
(described in the Fock space 3^) in the interaction picture with respect to the 
free dynamics of the fields |5ll7j . 
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4 Observables and instruments 



Let us consider now the case in which we are interested in the behaviour of 
the system Sj{, but any action on it is mediated by some input/output fields, 
for instance the electromagnetic field, represented in the Fock space 5". In this 
situation we can measure only field observables, from which we make inferences 
on S'k; we can speak of an indirect measurement. By choosing field observ- 
ables which commute also at different times in the Heisenberg picture, we can 
represent also measurements in continuous time. To this end we show how to 
construct such observables and how to eliminate the fields by a partial trace; 
in this way we obtain a description of the continuous measurement in terms of 
quantities (instruments) related to the system S:k alone |5I7) . 

To give the observables at all times, we have to give infinitely many com- 
muting selfadjoint operators or their joint spectral measure; the easiest way to 
do this is to work with the "Fourier transform" of such a spectral measure, the 
characteristic operator. The construction below involves the Weyl operators 

Hypothesis 3 (Elements of the characteristic operator). Let B^,B^,. .. , 
be commuting selfadjoint operators on Z and let us take c £ Ll^^(M.^;R^), 
b e L^„^(M+; Z) and h" € Lf^^{R+;Z), a = 1, . . . , m, such that 

lm{h°'{t)\hf^{t)) = 0, B°'h^{t)=0, Vt>0, Va,/3 = l,...,m. (16) 

Definition 6 (The characteristic operator). For any test function k = (ki, . . ., 
km) e L°"(M+;R'") let us define St(fc) e U{L^{R+;Z)), rt(fc) e L^{R+;Z) and 
the characteristic operator $t(fc) G 11(3^) by 

{St{k)f){s) = l(o,t)(s) [S(fc(s)) - a] f{s) + f{s) , V/ e L\R+; Z) , 

S{k{s))^l[e''-^-^^^\ ri(fc)(s) = l(o,,)(.s)r(fc;s), (17) 

m 

r(fc;s) =i^fc„(s)/i"(s)+ [S(fc(s)) -1] fe(s) , 



i ds Vfca(s)c"(s)+Im(6(s)|S(fc(s))6(s)) 

xW{rt{k);St{k)). (18) 

By dUl) and (HZ]), we get S(-fc) = S(fc)* and S(fc(i)) *r(fc; t) = -r{-k;t). 

Theorem 8 (0 Theorem 3.1]). Under Hypothesis\^ the characteristic operator 
introduced in Definitions^ has the following properties: 

1. localisation properties: $t(l(fj_42)fc) = (l(ti,t2)fc) ^ '^{'^{tiM))' ^ 
ti<t2< t; 
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2. group property: $t(0) =1, $t(fc)$t(fc') = $t(fc+fc'), Vfc, fc' € 

3. continuity: $f(Kfc) is strongly continuous in n £W and in t > 0; 
4- matrix elements: 

(e(5)|$.(fc)e(/)) = (e(g)|e(/)) 



X exp< / ds 



a/3 
m 

+ ij] &„(s)(c"(s) + (/."(s)|/(s)) + {g{s)\h-{s)}) 

a=l 

+ {g{s) + b{s)\{S{k{s))-TL) {f{s)+b{s))) 

5. $t(fc) is the unique unitary solution of the QSDE 

d$t(fc) = J2 Gy(t;fc)$t(fc)dA,,(t), $o(fc) = a, 
Goo(t;fc) =(fe(t)|(S(fc(t)) -1) 



a, 13=1 



Gjo{t;k) = (z,|r(A:;i)), Goj(i; fc) = -(r(fc; t)|S(fc(t))zj 
G.,(t;fc) = (z.|(S(fc(t))-l)z,). 



(19) 

(20a) 

(20b) 

(20c) 
(20d) 



Moreover, there exist a measurable space (f2, D), a projection valued meas- 
ure ^ on (51,1)), a family of real valued measurable functions |X(a,t; •), a = 
1, . . . , m, i > O} on il, a family of commuting and adapted selfadjoint operators 
{X{a, t), a — 1, . . . ,m, t > 0^ such that X{a, 0; uj) — 0, X{a, 0) = and, for 
any choice of n, ~ to < ti < ■ ■ ■ < tn < t, k^^ £ R, 

$t(fc) =expii^^4[X(a,tO-^(a,ii-i)] [ 

1=1 a=l ^ 

= / exp iV V«;^l(a,ti;L.)-l(a,ii_i;c^)] k(dc^), (21) 

^ l=lc.=l J 

where fc„(s) = ^27=1 («) '^L- 

From the unitarity and the group property we have $f(fc)* — $f(fc)~^ = 
^t{—k). Equation ((20)) is a right HP-equation with trivial initial space; it can 
be written also in left form, with the same coefficients. Note that Gij{t; ~k) = 
Gj^{t;k). 
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The observables X{a,t) can be identified by taking ka{s) — k, kp{s) = 
for /3 ^ a; then, the first equality in (PT|) gives e"''^^"'*) — $t(fc) and, by 
differentiation of the matrix elements p^ . we get 



e 



{g)\X{a,t)e{f)) = {e{9)\e{f)) / d,s{c"(.) 

+ {h%s)\f{s)) + (.g(s)|/i"(s)) + {g{s) + + h{s)])}. 



Let us choose the complete orthonormal system {z^, j = 1, . . . , d} in Z such that 
it diagonalises all the operators B^, . . . , i?™ and such that its first d! compon- 
ents, < d' < d, span the intersection of the null spaces of these operators; 
then, we have B°' = J2'i=d'+i-^i\^i)i^i\^ ^ write, on the 

exponential domain, 

d' ,t 



X{a,t)= f c«(s)ds + V / (hf{s)dA,{s) + h^{s)dAl{i 
Jo ,^1^0 ^ 

+ V / (dA,, + 6^dA,(s) + his) dAlis) + Ms)\^ ds) . (22) 

In quantum optical systems the continuous measurement of observables of the 
type J* (h"{s)dAi{s) + hf{s) dAl{s)j can be obtained by heterodync/homodyne 

detection, while terms like J* (^dAa + bi{s) dAi{s) + bi{s) dA|(s) + |fei(s)|^ ds^ 
are realised by direct detection, eventually after interference with a known signal 
if 6, ^ [31. 

A key point in the whole construction is that even in the Heisenberg picture 
the observables X{a,t) continue to be represented by commuting operators 
and, so, they can be jointly measured also at different times. Let U he a right 
unitary cocycle representing the system-field dynamics and define VT > the 
"output" characteristic operator by %^\k) := U{T)*^T{k)U{T). The key 
property giving the commutativity of the observables in the Heisenberg picture 
is $5,"* (l(o,i)fc) = $r*(fc): 0<t<T. This property follows from the fact that 
we have U{T) = U{T,t)U{t) (see Remark [1]), $t (l(o,t)fc) = ^t{k) and that 
U{T,t) e U{'K(8)Ji^t,T)) commutes with $t(fc) G U (?'(o,t)). 

Let s G §{3-C (8) 3^) be the initial system-field state. The characteristic func- 
tional of the process X (the "Fourier transform" of its probability law) is given 

by 

$t(fc) = Tr{$t(fc)C/(t)s(7(<)*} =Tr{$r'(fc)s}- (23) 

All the probabilities describing the continuous measurement of the observables 
X{a, t) are contained in $t(A;); let us give explicitly the construction of the joint 
probabilities for a finite number of increments. 

The measurable functions |x(a,t; •) , a — 1^ ... ,m, t > o|, introduced in 
Theorem [HI represent the output signal of the continuous measurement. Let us 
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denote by AX(ti,t2) = ta) - ^(1, ^i), • ■ • , ^(m, ta) - ^(m, ^i)) the vec- 

tor of the increments of the output in the time interval (ti, ^2) and by £,{dx; ti, ^2) 
the joint projection valued measure on R™ of the increments X{a, t2) — X{a,ti), 
a = 1, . . . , m. Note that, because of the properties of the characteristic oper- 
ator, not only the different components of an increment are commuting, but 
also increments related to different time intervals; this implies that the projec- 
tion valued measures related to different time intervals commute. Moreover, the 
localisation properties of the characteristic operator give 

^{A;ti,t2) = C{AX{ti,t2) eA)e t^)) , for any Borel set A C R'". 

(24) 

As in the last part of Theorem [U let us consider — to < ti < ■ ■ ■ < tn < t, 
ka{s) = J2?=i ^{t,-i.ti}{s) kL; then, we can write 



$t (fc) = Tr jexp i^i J2 [^("' - ^("' ) U{t)5U{tY | 



where the physical probabilities are given by 



,[AX(to,ii) e Ai,...,Al(t„_i,t„) e A„] 

= Trj (^f] ; ii-i, C/Wst/W* }■ 



Obviously, is the characteristic function of the physical probabilities 

[AX(<o, ii) e Ai, . . . , AX(i„_i, t„) e A„] and it uniquely determines them. 
The aim is now to reformulate the continuous measurement in terms of 
system alone, when the initial state is 

s = po®m))mf% PoeSm, feL^{R+;Z). (25) 

Let U (t) be a unitary, strongly continuous right cocycle and let us define 
U{t,s) by Eq. ([14]). Let $t(fc) be the characteristic operator introduced in 
Definition [6] under Hypothesis [3] and set 

$(fc;s,0 := $i(l(.,+oc)fc), 0<s<t. (26) 

By the definitions (|18l) . (|26l) and the points (l)-(3) of Theorem [5] one gets eas- 
ily $(0;s,t) = 1, $(fc;r, t) = $(fc; r, s) $(fc; s, i) and that $(fc;s,t) is strongly 
continuous in s and 

Definition 7. Let us take / G L^(M_|_; 2,) and < s < i. The reduced charac- 
teristic operator is the unique operator Gf{k; s,t) : £(5{) — )• £(J{) that satisfies, 

(v|^;/(fc;s,t)[X]?i) = (J7(t,s)«® V(/)|(^«)$(fc;s,t))C/(t,s)u(g) V(/)). (27) 
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Then, Tf{s^t) :— Qf(0]s,t) represents the reduced evolution operator for the 
observables of 

Theorem 9. In the hypotheses above, the family of linear maps Qf{k\s,t), 
t>s>0,fe i2(M+; Z), k e i°°(R+;R''), has the following properties: 

1. gf{k;s,s)=\- \\gf{k;sM < 1; 

2. Qf{k] s,t) is completely positive definite, i.e., for all integers n, test func- 
tions fc*, vectors (pi and operators Xi, one has 

n 

J2 {<t>^\Sf{k' - y;s,t)[X*X,]4,,) > 0; 

3. Qf{k;s,t) is a a-weakly continuous operator on £(IK) and it has a pre- 
adjoint Qf{k; s,t)^ acting on the trace class on Jf; 

4. for each X G L{J() the maps t 1— >■ Qf{k;s,t)[X], s 1— ?> Q f{k; s,t)[X] and 
K G /{Kk; s,t)[X] are continuous with respect to the a -weak topology of 

5. Vw,w e ■K, \JX e L{'K), 

{v\gf{k- s, t)[x]u) - {u{t, s, t))u{t, 

6. if f{x) = g(x) for all x £ {s,t), we get Qf{k;s,t) — Qg{k;s,t); then, 
Qf(k]s,t) is well defined for all f G L^q^(R; Z); 

7. gf{k;r,s)ogf{k;s,t)^gf{k;r,t), 0<r<s<t; 

8. for all s,t> we have Q f{k; s, s + t) — Q f^{ks;0,t) , where 

we have introduced the shifted functions fsix) = f{x + s), ks(x) ~ k{x-\-s), 
hs(x) — h{x + s), bs{x) — b{x + s), Cs{x) = c{x + s). 

Moreover, the evolution operator 'Tf{s,t), t > s > 0, f G L^(M+;2.), intro- 
duced in Definition^ enjoys the properties: 

(i) Tf{s,m]^l; r/(s,s) = l; \\rf{s,t)\\=l; 

(ii) Tf{s,t) is a a-weakly continuous operator on L{3{) and it has a pre-adjoint 

Tf{s,t)^ acting on the trace class on "K; 

(iii) Tf{s,t) is completely positive; 

(iv) for each X G £(!K) the maps t 1— > Tf{s,t)[X] and s i— > Tf{s,t)[X] are 

continuous with respect to the a-iueak topology of L{'K); 

(v) Tj{r, s) o Tf{s, t) ^Ts{r,t), < r < s < t; 
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(vi) if f{x) = g{x) for all x G (s,i), we have Tf{s,t) ~ Tg{s,t); then, Tf{s,t) 
is well defined for all f G ^^^^.(R; Z); 

(vii) for all s,t > we have 7/(s, s + t) = Tf^ (0, t), where fs{x) = f{x + s). 

Proof. The first statement of point (1) is immediate from the fact that ^{k; s, s) = 

1. The second statement follows from ||^^/(fc; s, t)[X]|| < X®$(A;;s,<) = 

\\X\\; the first step is from the definition (P7|) . the unitarity of U{t,s) and the 
normalisation of the coherent vector ip{f), while the second step is due to the 
unitarity of the characteristic operator. 

By using $(A;* — s) — $(— s)* $(— F;*, s) and the definition of 
Qf{k]S,t), one gets immediately 



^ {^^\gf{¥-k=;s,t)[X*X,^,) 



.$(-/c-';t,s)C/(i,s)0,(8)V(/) 



> 0, 



which is point (2). 

Any T e 'J(J{) can be written as r = \un){vn\ for some choice of the 
vectors u„, w„ in 3-C. Then, we have 

Trjc {Gfik; s, t)[X]T} = ^(«„|e/(fc; s, t)[X]ur,) 

n 

= J2 {U{t, s) vn ® I {x ® $(fc; s, t)) U{t, S) Un ® Hf)) 

n 

{ [X ® $(fc; s, t)) U{t, s) (r (g> (^(/)|) U{t, s)*} 

^:TT^{Xgf{k;s,tU[T]}, 



Tr 



which defines the pre-adjoint. The existence of the pre-adjoint of Qf{k]s,t) 
implies its a-weak continuity [231 Corollary of Theorem 1.13.2, p. 29] and this 
completes the proof of point (3) 

By point (1) Qf{K,k] s, t) is bounded uniformly in s, t and n. By Proposition 
1.15.2 in 24 , the weak and the cr-weak topologies are equivalent on the bounded 
spheres; so, it is enough to prove the weak continuity. Let us set 01 := s) v<S) 
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^(/), := U{t,s)u(E)ip{f), X := X (E)<^{k;s,t). Then, we have 
livlGfik; s,t + e)[X]u) - {v\gf{k; s,t)[X]u)\ 

= {U{t + e,t)(f>i\X^{k;t,t + e)U{t + e,t)c^2) - (^il^^a) 



< 



((C/(i + e, t) - 1) </>! |X$(fc; t, i + e)J7(t + e, t)02> 



61 i + e) (C/(t + e, i) - 1) 02> 



L|x($(fc;i,t + e)-l 



< 



X$(fc; t, t + e)J7(t + e, ||</.2|| || {V{t + £,<)- I) 0i|| 



+ 



+ e) 



||0i||||([/(< + e,O-l)02|| 
|l0i|| II + e) - U) 02|| < ll^ll { \\A\ W(t + e,t)- I) 

+ ||«|| [||(c/(t + e,t)~ 1) 02II + II ($(fc; + e) - a) ^2]] }, 



which gives the continuity in t. The continuity in s can be proved in a similar 
way. By similar steps we get 

\{v\gfiK'k;s,t)[X]u) - {v\gfiKk;s,t)[X]u)\ 

< \\X\\ \\v\\ f s, t) - $(Kfc; s, t)) U{t, s)u ® 



which gives the continuity in n, due to point (3) in Theorem [8l This ends 
the proof of point (4). Points (5) and (6) are immediate by the localisation 
properties. 

By using the identification = V'(/(o,r)) ® '4>{f{r,s)) ® i^U(s,t)) V'C/c*) 

and the localisation properties $(fc;a, 6) € lt(3^(a.b)), U{b,a) S (g) 3^(0,6)), 
we have 



{v\gf{k;r,s)ogf{k;s,t)[X]u) 
= (;7(s,r) I {gfik;s,t)[X](g>^{k;r,s)^ U{s,r) (u ® 

= s) [U{s, r) {v (g> iPifir^s))) ® I 

((X ® $(fc; s, i)) » *(fc; r, s)) {/ (t, s) [U{s, r) {u ® ® ^{hs,t))\ ) 

= (U{t, r) [v ® I ® ^' ^)) ® ^(/(r.t))) ) 

= (Hg/(fc;r,t)[X]u), 
which gives point (7). Finally, by Eqs. (P, (HI]), ([15]), (HH) we have 

{v\gf{k-s,s + t)[X]u) 

= (C/(i)v ® I ® e,$(fc; s, s + t)Ql) U{t)u ® 



{U{t,Q)v®i}{fs)\ [X®^ks-M 



[/(t,0)u(»V(/s)), 
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and point (8) follows. 

By the particularising the previous statements to the case fc = 0, we get the 
properties of the evolution operator. □ 

The definition of the reduced characteristic operator has been given in such 
a way that it is sufficient to construct the characteristic functional (1^51) when 
the initial state is given by Eq. ([25]): $t(fc) = Tr {^/(fc; 0, t)[ll]po}- So, the 
reduced characteristic operator determines all the probabilities of the output. 
However, the reduced characteristic operator gives something more: the states 
after the measurement, conditional on the observed output. This is obtained 
through the correspondence with the instruments representing the continuous 
measurement, see pp. 244-245] and [5]. 



5 The evolution equations 

Up to now, we have only made use of the cocycle properties of J7(t), but we 
are interested in finding the infinitesimal generator and the evolution equation 
of the reduced characteristic operator and for that we need also the QSDE for 
U (t) . The reduced characteristic operator comes out from the product of three 
terms: the operators $f(fc), U{t) and U{t)*. To compute the differential of this 
product we have to use two times the second fundamental formula of quantum 
stochastic calculus. 

Our first step will be to differentiate the unitary process 

*t(fc) := J7(i), fee (28) 

then, we shall use the second fundamental formula of quantum stochastic cal- 
culus to elaborate the expression giving the reduced characteristic operator. 

Lemma 10. Let Hypotheses[Il\^\^hold and the functions c{t), b{t), h°'{t) be 
locally bounded in time. Then, ^'t(fc), defined by ([28]), can be expressed as the 
quantum stochastic integral on Dom(C^/^) £ 

*t(fc) = ]l+^ / (1® $,(fc))Af,,(s;fc)[/(s)dA,y (s), (29) 

where 

d . 

Mooit; k)^K + ^(r(-fc; t)\zr)Rr + <^ {b{t)\ {S{k{t)) - l) b{t)) 

+ i^fc„(i)c"(t)-- K{t){h^{t)\h^{t))kp{t)W, (30a) 

a=l a,/3 = l ^ 

d 

Moj{t; k) = iV, + ^(r(-fc; t)\z,)Srj , j > 1 , (30b) 

r=l 
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MM{t;k)^Y,{''4^i'^(''^))''r)Rr + {z;\r{k;t))l, i>l, (30c) 

r>l 

M,,{t; k)^J2 (^^|S(fc(^))z,,) Srj - S,jt, I, J > 1, (30d) 

with Dom(Afoo(t;fc)) = Dom(/-i:), Dom(Afoj (i; fc)) = I)om{Nj) D Dom(i^), 
Dom(M,o(t;fc)) D flLi Dom(i?fc) D Dom(is:) U Dom(i^*), Doni(A%(t; fc)) = 
i,j = 1, ■.■,(!. 

Moreover, yf,g G M and\/u,v G Dom(C^/'^), one has 
{U{t)v®e(g)\ (X®^t{k)) U{t)u (g> e{f)) = {v\Xu){e{g)\e{f)) 

+ V / ds^|(C/(s)i;®e(g)|(X®$,(fc))A/,,(s;fc)C/(s)M®e(/)) 

+ {Fj,U{s)v (g) e(5)|(X ® $,(fc))C/(s)M® e(/)) 

+ {Fi^Uis)v (g> e{g)\{X (g> $,(fc))M,,(s; k)U{s)u (g> e(/))}/,(s). (31) 
i>i 

Let us recall the convention fo{s) = .go(s) = 1- 

Proof. By Eqs. ([8]) and (|20t , the second fundamental formula of quantum stochastic 
calculus, ^t(-k) = ^t{k)* and Gji{s;-k) = Gij{t;k), we get for /,g S M and 
u,v€ Dom(Ci/2) 

(z; ® e{g)\^t{k)u (g> e(/)) - («|«)(e(5)|e(/)) 

= J2 f ds'^{v(g>iti-k)e{g)\M,,{s;k)U{s)u(^eif))f,{s), (32) 

where 

Afy (s; /c) := i^,^- + Gy (s; fc)! + ^ G,;^ (s; k)Frj . (33) 

By inserting the explicit expressions of the elements of the matrices F and G 
into Eq. ([33|) we get Eqs. (|30|) . The statements about the domains follow from 
Hypothesis [T] point (iv), Proposition [T] point (1) and the fact that the operators 
Sij are bounded. 

It is easy to check that the processes (ll(8)<&s(fc)) Af^^ (s; k)U (s) are stochastic- 
ally integrable, by using the fact that ^s{k) is unitary, the functions Gij{s;k) 
are locally bounded, due to the boundedness assumption on c, b, h", and the 
processes FijU(s) are stochastically integrable by hypothesis. Then, Eq. p9p 
follows from Eq. (|32|) and the first fundamental formula of quantum stochastic 
calculus. 

By the second fundamental formula of quantum stochastic calculus applied 
to {X* (g) !)[/(<) and ^'((A:) we get immediately Eq. □ 
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For A, r G Z (with components denoted by Xj and rj) let us define the 
operators 

d 

B,{\) i?, + ^ S,j\j, i = 1, . . . , d, (34a) 

d II , ||2 d 

ij—l i=l 

By taking into account Hypothesis [1] and Proposition [1] we have 

Dom(B,(A)) = Dom(i?,) D Doni(iv:) U Doni(iv:*), 
Dom(/i:(A, r)) = Dom(X) D Dom(C^/2). 

Again by Hypothesis [1] and Proposition [U the domains of the adjoint of the 
previous operators contain Dom(i^*) D D and on Dom(i^*) we have 

d 

B,{X)* =R*+Y,^S,*, z = (35a) 

d II . ||2 d 

K{\ ry^K*-J2 '>^S*^R, - ^ 1 + ^ r, B,(A)*. (35b) 
Finally, for a^, c G R"', b e Z, h e Z"^ we define also 

m I m 

C(k,6,c,/i) := (6|(S(K)-l)6)+i^K„c"-- ^ Kc.(/i"|/i'^)k/3 . (36) 

a=l a,/3=l 

Proposition 11. Let Hypotheses [TJ [H hold and the functions c{t), b{t), 
h^'it) be locally bounded in time. Then, V/ G M, Vfc e L°°(R+;M™), Vit, w G 
Dom(C^/^) we have 

{U{t)v ® ^(/)| ® $t(fc)) C/(t)u ® V(/)) = 

+ ^ ds |({/(s)« ® ^{f)\{X ® $s{k))K{f{s), r{-k, s))U{s)u <E> V(/)) 
+ {K{f{s), r{k, s))U{s)v (g, ij{f)\{X (g> ^s{k))U{s)u (g, ij{f)) 

d 

+ ^(z,|S(fc(s))z,)(B,(/(s))[/(s)z;®V(/)|(^®$.(fc))i?,(/(.'^))f/(s)w®^(/)) 

+ C{k{s),b{s),c{s),h{s)){U{s)v ® V(/)| ® ^s{k))U{s)u ® ^(/))|. (37) 

Proof. The statement fohows by direct computations, by inserting the explicit 
expressions of , Mij{t; k), Nj into Eq. (PT|) with g ^ f ■ □ 
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Let 2) C be the linear span of the rank-one operators of the type 

|'!/')(<?!'| with tp, (j) G Dom(F*). By using the operators (l35t . we define, Vip, cj) G 
Dom(F*), Vu, V e^K, 

{vllCfitmmu) = {v\i,){K{f{t), r{~k,t)ycj,\u) + {v\K{f{t), r(fc,t))»(0|u) 

d 

+ {z,\S{k{t))z,){v\B,{f{t)y,/j){B,{f{t)yct>\u} 

+ C{k{t),b{t),c{t),h{t)){v\i^){<l>\u); (38) 
then, by hnearity, we extend /C^(<) to D. 

Corollary 12. Let Hypotheses[J[[^[^ hold and the functions c{t), b{t), h°'{t) 
be locally bounded m time. Then, V/ G M, Vfc G i°°(]R+; M"), Vm, u G 'K, 
yX G CD, we have 

{v\gf{k;0,t)[X]u)^ {v\Xu)+ [ {v\gf{k;{),s)[lC){s)[X]\u)ds. (39) 

Jo 

Proof. By using the notations above, Proposition[TT] gives immediately Eq. (|39p 
Vu,D G Dom(C^/^). Being X G D, the operator /Cy(s)[X]] turns out to be 
bounded; moreover, we have ||t?/(fc; s, t)|| < 1. Then, we can extend (|39| to any 
u, w G J{. □ 

By introducing the pre-adjoint of Gf{k\ 0, t) and extending to the whole 
trace class we get: VX G CD, Vt G T(5f), 

TT^{Xgf{k-Q,tUT]}^Tr^{XT}+ [ Tr.^{/Cj(s)[X]g/(fc;0,s)4r]}ds, 

"'0 

(40) 

with initial condition Qf{k;0,Q)^ ==1. For A: = and r G S(J{), Eq. pi)) 
is a quantum master equation and the formal pre-adjoint of IC'j{t) is known as 
Liouville operator. We can say that ()40p is a generalisation of a quantum master 
equation, which includes the continuous measurement. 

The problem which remains open is to prove the uniqueness of the solution 
of Eq. ([55)1 or of Eq. We note that in the case of quantum dynamical 

semigroups the positivity plays a role in the analogous problem, see Theorem 
m while in the case of Eq. ((39)) we have only that Q is positive definite in k. 

6 An example: the degenerate parametric oscil- 
lator 

The degenerate parametric oscillator is the physical system which was used to 
produce squeezed light. |19l26l27j The squeezing of the light was revealed by bal- 
anced homodyne detection, a measurement scheme which is indeed described 
by continuous measurements of diffusive type. [416] . Such a quantum optical 
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system is constituted by an optical cavity closed by two partially transparent 
mirrors with inside a crystal with a x'^^'' non-linearity. Only two cavity modes 
of the electromagnetic field inside the cavity are relevant: the subharmonic field 
of frequency uic (a quantum oscillator with annihilation and creation operators 
a, a^) and the pump field of double frequency (with annihilation and creation 
operators 6, &^). The pump field is populated by a resonant laser entering the 
cavity, the crystal couples the two modes and the light coming out of the cavity 
is detected by homodyne devices and/or photocounters. The degenerate para- 
metric oscillator is well studied from the point of view of theoretical physics 
and quantum optics in [11] Chapts. 9, 10, 12, 18]. Here we want to prove that 
the mathematical model of the degenerate parametric oscillator with direct and 
homodyne detection can be rigourously formulated and gives an example of the 
theory we have developed. 

The formal master equation is given in [IH Eq. (9.97)] and reads 

p{t) = -i[Ho,p{t)] - i [Ae-^'"^*6t + Ae2-^*&, p{t)] 

+ K (n + 1) {2ap{t)a'' - a^ap{t) - p{t)a^ a) 

+ Kn (2aV(0a - aa'^ p{t) - p{t)aa^) + KpUp {2b'' p{t)b - bb'' p{t) - p{t)bb'') 

+ Kp {rip + 1) {2bp{t)b^ - b''bp{t) - pitpb) . (41) 

The Hamiltonian term Hq contains the free energies of the modes and the in- 
teraction due to the x^^) non-linearity: 

Ha = wca^a + 2ujcb^b + |^ (a^^^ - b^a^) . (42) 

For the constants we have ujc > 0, 5 G K, g ^ 0, k > 0, n > 0, Kp > 0, 
rip > 0. 

This model, plus detection, can be rigourously formulated in the set up 
developed before. First, the Hilbert space is identified with the span of the 
eigenvectors of the two number operators and the creation and annihilation 
operators are defined. Let us take the Hilbert space 5f = £^(N) (8)£^(N) with 
canonical orthonormal basis {en,m, n,m> 0}. The creation, annihilation and 
number operators for the subharmonic mode are defined by 

Dom(a) = Dom (a^) = < u G !K : n < +00 > , (43a) 

[ n,m>0 J 

en,m = Vn + 1 Cn+i.m, a eo,m = 0, ae„,m = V"e„-.i.m, if n > 0, 

(43b) 

Dom (a^a) = < u G J{ : n"^ \un,m\^ < +00 > , a) aen,m = nen^m- 

I n,m>0 J 

(43c) 

An analogous definition holds for the operators b\ b, b''b, which act on the 
second factor of the tensor product. 
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In constructing the model we have to reproduce the effects contained in 
the master equation (|41l) and to introduce the measurement. So, we have to 
introduce losses at the mirrors and thermal dissipation in the crystal and at the 
walls of the cavity. We have also to introduce the possibility of injecting laser 
light feeding the pump mode. Moreover, we consider the direct detection of 
photons with two photocounters, chosen one to be sensible to photons around 
frequency ojc and the other to photons around frequency 2ujc- Finally, we 
consider homodyning around frequency ujc- To realise all these features in the 
mathematical model we need many channels, but some channels with similar 
structure can be collected together and the minimal number is d 8. We 
use channels 1 and 2 to describe the light reaching the two photocounters and 
channel 3 for the light reaching the homodyne detector, channel 4 is the one 
used for the injection of the laser, channels 5-8 describe losses and thermal 
dissipation. There is no scattering effect which mixes the channels. The channel 
operators and the unitary matrix of system operators we need are 

i?i = Pi b, R2 = ai a, i?3 = a2 a, R4 = /32 b, R5 = /Jg 6, (44a) 

Rq — as a, i?7 — /34 b\ Rs — Q!4 a^, Sij — (5y 1, (44b) 

|ai|^ + |a2|V|a3|^ =2K(rT+l), {aif ^ 2Kn, (44c) 

l/3i|' + |/32|' + |/?3|' =2«;p (rip + 1), I/34I' = 2Kpnp. (44d) 

The operator K has to include the Hamiltonian Hq and to satisfy Eq. ([TU)) : 
so, it must have the formal expression 

1 ^ 

K = —iHo — - Ri* Ri = '- (^a^^b — 6^a^) — (ktI + Hpfip) 1 

i=l 

- (iwc + K (2n + 1)) a) a - (2iwc + Kp (2np + 1)) 5^6. 
Rigourously, by defining u„^,„ = if n < and/or m < 0, we have 

Ku = y2ku{n,m)en,m , Dom(A') = < u : |fc„(r7,, m)|^ < +00 > , (45a) 



K{n,m) ^ n(n - l)(m + 1) u„-2,m+i - \J m{n + l){n + 2) u„+2,m-i 

— [ku + KpHp + icjpn + K {2n + l)n + 2ia;c'm + Hp (2np + 1) m] Un,m ■ (45b) 

Theorem 13. Let us construct the F -matrix by setting Fqq = K, F^q — Ri, 
Fqj = Nj =: —Rj, Fij = 0, i,j > 1. Then, Hypothesis]^ hold true with D = D 
given by the linear span of the basis {en,m, n,m>0} and with C = , where 
TV flta + 26^6. 

Proof. By applying the definition of adjoint and Riesz lemma j22| one can easily 
check that a* = a\ a'^* = a and the same for 6, b\ as it is well known. In 
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particular all operators Ri, Ni are closed [13] ■ By (I43ap we have the domain 



Db.n := Pi Doni {F^j) ^ I u "K : ^ {n + m) |u„,„|^ < 



-oo 



i, j I n.m>0 

(ij)#(0.0) 

(46) 

Again by the definition of adjoint, [22^ we get K*, which turns out to be 
defined by Eqs. P5|) with the substitutions ujc —^c^ 9 ~^ ~9- From the 
definitions of K and K* we get, by direct computations, the dissipativity con- 
ditions (|10l) and K** = K. In particular also K is closed. 

For every u ^ D we get, from the dissipativity condition 



\Kuf > (Ku 



u \ / u 



_ \{Ku\u)\' ^ iRe{Ku\u)) 



2 (j:,\\RkU 



uf - \\uf A\\uf 



which gives ||i?feu||* < 4 and Dun D Dom(_ft'). Analogously, we get 

Djipf D Doni(ii'*). Up to now we have proved conditions (i), (ii), (iv), (v), (vii). 

To prove condition (iii) we have to show that the set D given in the Pro- 
position is a core for K, K* , a, a^, b, 6^. By {u\aen,m) = \pn Wn-i,™, we 
get 

{m S J{ : 30 S ;H : (u|ae„.„i) ~ (0|e„.„i) Vn,rn} = Dom(at). 

Similar considerations hold for a^, 6, 6^ and I? is a core for a, a^, 6, fe^. Ana- 
logously, by writing the expressions of Ken,m e K*en,m we deduce that D is a 
core for K and K* . 

By (Uni) both K e K* are dissipative P^, Definition 4.1 p. 13] and, by [H 
Corollary 4.4 p. 15], they are generators of contraction semigroups in IH. This 
complete the proof of point (vi). 

Finally, let us consider point (viii). 

We have Dom(Ci/2) = |m € J{ : E„,™>o(" + 2m)'* < +oo|, which 

is obviously contained in Dom(F) = Dom{K) given in (pS)) . 

^ 1/2 ~ 

For any e > take = D = D. Then, Ce D^<Z D because is diagonal 
in the canonical basis. For large n and/or m, en,m goes as l/(e(?^ + 2m)^) 

1 /2 

and each operator F*^Ct \d^ is bounded (the worst case is for _F'q*q = K*). This 
is point (a). 

By explicitly computing the left hand side of the inequality in point (b), we 
see that we have to prove the inequality 

2Re^(u, + ii?,uo|[Ce,i?.]uo) < Yl {bl{'^^\C.u^) + 62 llu.f) . 

The proof of this inequality is long and we give only a sketch. 

Let us note that for any function of the operator N one has af{N) = f{N + 
l)a, fiN)a\ 6/(iV) = /(TV + 2)6, /{N + 2) ^ f{N)b^ . We 

take fix) = (YT^ ^^^^ = ^ 2; < 0, so that f{N) = C^. 
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By these relations and standard estimates one can obtain 

Re{aau\[C„aa]u) < \af {a^au\ \f{N-l) - fiN)\u). 

By the relations above, expansion in the canonical basis and standard estimates 
we can prove also 

2 Re{v\[C,,aa]u) < {a^au\ |/(7V - 1) - /(iV)| u) + {v\ \f{N) ~ f{N + l)\v). 

Analogous estimates can be obtained in the cases involving a\ b, . All together 
these results give 

2ReJ2{u^ + I R^uo\[C,, R^]uo) < J2 {u.\\f{N + 2)-f{N)\u,} 

i>l 1=1,4,5 

+ Ak (n+l) {a^auo\ |/(iV - 1) - f{N)\uo) + {a^aus\ \f{N - 1) - fiN)\us) 
+ 4Kn{{a^a+l)uo\\f{N + l)-f{N)\uo)+ {u,\\f{N + l)-f{N)\u,) 

■i=2,3,6 

+ Akj, {np + 1) {b^buo\ \f{N - 2) - f{N)\ uo) + {b^bm\ \f{N - 2) - /(iV)| uy) 

( {b^b + l) uo| |/(^^ + 2) - f{N)\uo). 

By using the specific form of / and 2{x ~ 1) > x for x > 2 and 3 (a; — 2) > x 
for a; > 3, we get 



1 



fix + l) 



< 



15 



X 



1 



fix + 2) 



64 

< — , for a; > 2; 

X 



/(2:-2) 



< 



80 



< 



for X > 1 
648 



for a; > 3. 



These inequalities can slightly modified to include also the cases x — 0,1,2. 
Then, by further straightforward estimates, one gets 



2RcJ2{u^ + \ R^u^C,,R,\u^ < 4 (Kn + 16Kpnp) ||Mof + 324(u7|CeM7> 
+ 64(m8|C,?/8) + (||u4' + 15(?/.|CeU,)) +16 Y (h4' + 5(u,|C7eU,) 

1=2,3,6 1=1,4,5 

+ 8 (32K(n + 1) + 162Kp(rip + 1) + 15Kn + 60Kpnp) (uo|Ceito). 
This ends the proof of the inequality. □ 

To use the number operator N — a^a + 26^6, which commutes with Hq, is 
suggested by [101 Chapter 3] , where the conservativity property of the minimal 
quantum dynamical semigroup in a similar model is proved. 

Proposition 14. Hypothesis\^ holds for the model of this section. 
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Proof. We prove the sufficient condition of Tfieorem [S] with A = K* , Lk = 
= —Rk- Let D be as in Theorem 1131 and, on their maximal domains, let us 
introduce the operators 

Q :~ wN + 2Kn + 2Kpnp , w := max {2k {2n + 1) , Hp (27Tp + 1)} , 

Z :— 2k (2n + 1) a^a + 2Kp {2rip + 1) 6^6 + 2Kn + 2Kpnp . 
By defining also v := min {2k {2n + 1) , Kp {2np + 1)}, on Z? we have 

V w — V 

< — Q H {2Kn + 2Kpnp) < Z <Q. 

w w 

In particular we get Dom(Q) — Dom{Z) ~ Dom(A^) — Dom(a^a) n Dom(6^fe), 
D C Dom(Q) C Dom (Q^/^) = Dom (Z^^'^)- The set is a core for Q^/^ 

In the proof of Theorem [T3l it is shown that Dom{K*) C Db,n- But one 
can check that Drn = Dom((5^/^), so, we have Dom(/^*) C Dom (Q^/^) = 
Dom(Zi/2). 

Finally, we get nfe>i -Rfc(Dom(iir*^)) by the fact that the i?fcS are propor- 
tional to a, a\h or &^ and that T)ouy{K*'^) C Dom(ata) n Dom(6t), as one can 
check. 

For u € D we get by direct computations 



-2Re{u\K*u) = \\Rk 



l|2 
u\ = 



k>l 



z'/\ 



z'^\ 



< 



w\\au\\ +2w\\bu\\ + 2 {ku + KHp) \\u\\ 



2Re{Q^/^u\Q^/^K*u} + ^ \\Q^/^RkU 



k>l 



2 (ku + KHp) \\u\\ —2wK\\au\\ —4:WKp\\bu\\ 



< 



Then, these inequalities can be extended to the domains required in Theorem [5] 
and this ends the proof. □ 

In order to describe the two photocounters and the homodyne detector we 
have to specialise the observables ((22|) : what we need is to take m — 3 and [5] 



X{a,t) 



-^aa (^) ; 



a = 1, 2, 



This means that the quantities in Hypothesis [3] are 

c(0 = 0, 6^=0, Bi = |zi)(zi|, B^ = \z2){z2\, B^=0, (47a) 
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h\t) = h''{t) ^ 0, hl{t) = 5,3e'(«^-"^*). (47b) 

This choice trivially satisfies Hypothesis [3] and the expressions of the quantities 
in Definition ini become 

2 

S(fc(s)) =l + ^(e'^^(^)-l) \z,){z,\, rt{k){s)^l^o,t){^)ih{s)h^{s), 

3 = 1 

2 

(St(%)(s) - (e"=^(^) - l) g,(.)z, r(fc; s) = ifc3(s)/i=^(s). 

Finally, in order to describe a coherent monochromatic laser pumping the 
6-mode as in the source term in the master equation (j4ip . we have to take a 
coherent state of the field with /-function given by 

iAe~2''^c^* 

m)^5,A - l(o.T)ft)- (48) 

P2 

We are assuming /32 7^ and we understand that T is a large time (needed to 
have an L^-function), but that T — > +oo in the reduced characteristic operator. 

In conclusion the model just described is well defined, as it satisfies all the 
hypotheses introduced in this paper. Moreover, one can check that the asso- 
ciated formal master equation (Eq. (|40|) for k — Q) reduces to Eq. (|4T|) . as we 
wanted. 
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